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Abstract 

Liouville conformal field theory is considered with conformal boundary. There is a fam- 
ily of conformal boundary conditions parameterized by the boundary cosmological constant, 
so that observables depend on the dimensional ratios of boundary and bulk cosmological 
constants. The disk geometry is considered. We present an explicit expression for the expec- 
tation value of a bulk operator inside the disk and for the two-point function of boundary 
operators. We comment also on the properties of the degenrate boundary operators. Pos- 
sible applications and further developments are discussed. In particular, we present exact 
expectation values of the boundary operators in the boundary sin-Gordon model. 

1. Liouville field theory 

During last 20 years the Liouvlle field theory permanently attracts much attention mainly 
due to its relevance in the quantization of strings in non-critical space-time dimensions [Q 
(see also refs.[0, ^ Q). It is also applied as a field theory of the 2D quantum gravity. E.g., the 
results of the Liouville field theory (LET) approach can be compared with the calculations 
in the matrix models of two-dimensional gravity |^, ^ and this comparison shows 0, ^ that 
when the LET central charge c^, > 25 this field theory describes the same continuous gravity 
as was found in the critical region of the matrix models. Although there are still no known 
applications of LET with < 25, the theory is interesting on its own footing as an example 
of non-rational 2D conformal field theory. 
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In the bulk the Liouville field theory is defined by the Lagrangian density 

C = ^ida<PY + l^e'"^ (1.1) 

where (f) is the two-dimenesional scalar field, b is the dimensionless Liouville coupling constant 
and the scale parameter fi is called the cosmo logical constant. This expression implies a 
trivial background metric Qab = Sab- In more general background the action reads 

^buik = ^ / b^'^aM^ + QR<P + A-Kfie'''^] ^d'x (1.2) 
r 

Here R is the scalar curvature associated with the background metric g while Q is an im- 
portant quantity in the Liouville field theory called the background charge 

Q = b+l/b (1.3) 

It determines in particular the central charge of the theory 

CL = l + QQ^ (1.4) 

In what follows we always will consider only the simplest topologies like sphere or disk 
which can be described by a trivial background. For example, a sphere can be represented 
as a fiat projective plane where the fiat Liouville lagrangian ( |1.1| ) is valid if we put away all 
the curvature to the spacial infinity where it is seen as a special boundary condition on the 
Liouville field 

(j){z,z) = -Q\og{zz) + 0(1) at \z\ ^ oo (1.5) 

called the background charge at infinity. 

The basic objects of LFT are the exponential fields Va{x) = exp(2a(/)(x)) which are 
conformal primaries w.r.t. the stress tensor 

T{z) = -{d<pf + Qd^(j) (1.6) 
f{z) = -id(t>f + Qd^(l> 

The field Va has the dimension 

A, = a(Q-«) (1.7) 

In fact not all of these operators are independent. One has to identify the operators Va and 
Vq_q, so that the whole set of local LFT fields is obtained by the "folding" of the complex 
a-plane w.r.t. this refiection. The only exception is the line a = Q/2 + iP with P real where 
these exponential fields, if interpreted in terms of quantum gravity, seem not to correspond 
to local operators. E.g. in the classical theory they appear as hyperbolic solutions to the 



Lioville equation and "create holes" in the surface R ITO[. Instead, these values of a are 
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attributed to the normalizable states. The LFT space of states A consists of all conformal 
families [vp] corresponding to the primary states \vp) with real < P < oo, i.e., 



^=09 H (1.8) 

PG[0,oo) 

The primary states \vp) are related to the values a = Q/2+iP and have dimensions 

while other values of a are mapped onto non-normalizable states. This is a peculiarity of the 

operator-state correspondence of the Liouville field theory which differs it from conventional 

CFT with discret spectra of dimensions but make it similar to some conformal cr-models 

with non-compact target spaces. In what follows the primary physical states are normalized 

as 

{vp,\vp) = 7i6{P - P') (1.9) 

The solution of the spherical LFT amounts to constructing all multipoint correlation 
functions of these fields, 

{Xi,... ,Xn) = {VaiiXi) . . .Vaj^{xN)) (1-10) 

In principle these quantities are completely determined by the structure of the operator prod- 
uct expansion (OPE) algebra of the exponential operators, i.e. can be completely restored 
from the two-point function 

{Va{x)VM) = (1-11) 

which determines the normalization of the basic operators and the three-point function 

Gai,a2,aaiXl,X2,X3) = Ixial^"^' |a;23|^'^' |X31 j^'^' C(ai, ^2, Q^s) (1-12) 

with 7i = Aa^ - - A„3, 72 = - Aa.^ - A^-,, 73 = - A^^ - Aq2 . Once 
these quantities are known, the multipoint functions can be in principle reconstructed by 
the purely "kinematic" calculations relied on the conformal symmetry only. Although these 
calculations present a separate rather complicated technical problem, conceptually one can 
say that a CFT (on a sphere) is constructed if these basic objects are found. 

For LFT these quantities were first obtained by Dorn and Otto|]ll], |12| in 1992 (see also 



||13|). We will present here the derivation of the simplest of them, the two-point function 
D{a), to illustrate a different approach to this problem proposed more recently by J.Techner 



Ij] which seems more effecticient. Close ideas are also developed in the studies of LFT by 
Gervais and collaborators [^. We will use similar approach shortly in the discussion of the 
boundary Liouville problem. 

Among the exponential operators Va there is a series of fields V-nb/2, n = 0,1, . . . which 
are degenerate w.r.t. the conformal symmetry algebra and therefore satisfy certain linear 
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differential equations. For example, the first non-trivial operator V-b/2 satisfies the following 
second order equation 



^d^+T{zyjV^b/2 = (1.13) 

and the same with the complex conjugate differentiation in z and T{z) instead of T. In 
the classical limit of LFT the existence of this degenerate operator can be traced back to 
the well known relation between the ordinary second-order linear differential equation and 



the classical partial-derivative Liouville equation [jT6|. The next operator V^h satisfies two 



complex conjugate third-order differential equations 

(J^d^ + 2Tiz)d + (1 + 2b^)dTiz)^ V.b = (1.14) 

and so on. It follows from these equations that the operator product expansion of these 
degenerate operators with any primary field, in the present case with our basic exponential 
fields V^, is of very special form and contains in the r.h.s only finite number of primary fields. 
For example for the first one there are only two representations 



V-l,/2V^ = C+ [V„_fc/2] + C_ [K.+,/2] (1.15) 

where C± are the special structure constants. What is important to remark about these 
special structure constants is that the general CFT and Coulomb gas experience suggests 
that they can be considered as "perturbative" , i.e. are obtained as certain Coulomb gas (or 
"screening" ) integrals |T7| , . For example in our case in the first term of ( p..l5| ) there is no 



need of screening insertion and therefore one can set C+ = 1. The second term requires a 
first order insertion of the Liouville interaction —fi J exp[2b(j))d'^x and 



C- = f/'x(r,(0)\/_,/2(l)e2^*(^VQ_„_V2(oo)> 

= ^{-b^H2ba) ^'-''^ 

where as usual 7(x) = r(a;)/r(l— x). It is remarkable that all the special structure constants 
entering the special truncated OPE's with the degenerate fields can be obtained in this way. 

Now let us take the two-point function D{a) and consider the auxiliary three-point 
function 



( 14 (Xi ) Va+b/2 {X2 )V^b/2{z)) 

Then, tending z — xi we see that in the OPE only the second term survives and in fact our 
auxiliary function is ~ C-D{a + b/2). Instead tending 2; — >■ X2 we can "lower" the parameter 
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of the second operator down to a which results in Cj^D[ol). Equating these two things we 
arrive at the functional equation for the two-point function 

F^T^ = C_ a 1.17 

This equation can be easily solved in terms of gamma-functions 



D{a) = (7r/i7(6^)) 



{Q-2a)/b 



7 (26a - 



6^7(2 -2a/6+ 1; 



;i.l8) 



which coincides precisely with what was obtained for this quantity in the original studies. 

In fact there are many solutions to the above functional equation. It is relevant for the 
moment to stop at the remarkable duality property of LFT. Besides the abovementioned 
series of degenerate operators V-nb/2 there is a "dual" series with b replaced by 1/6. This 
results in another "dual" functional equation for D{a) with the shift by 1/6 instead of 6. 
The solution becomes unique (at least if these two shifts are uncomparable) Jill. Note that 
these two equations are compatible only if in the dual equation the cosmological constant fi 
is replaced by the "dual cosmological constant" /2 related to fi as follows 



7r/i7(l/62) = (7r/i7(62)) 



;i.i9) 



With this definition of fi the duality property, which turn out to hold exactly in LFT, can be 
formulated as the symmetry of all observables w.r.t. the substitution 6 — 1/6 and /i — /t. 
The same way one can readily obtain and solve the functional equations for the three- 



point function [14] which reads 



C(ai, q;2, as) = 71/27(6)6 



(,2M,2-2fe2 



;i.20) 



ToT(2ai)T(2a2)T(2a3) 



T(ai + a2 + as — Q)T{ai + a2 — a3)T(a2 + as — ai)T(as + ai — a2) 
where a special function T(x) has to be introduced 



log T(x) 



°° dt 
T 



sinh^ {Q/2-x)t 
sinh(6t) sinh(t/6) 



'1.211 



This integral representation is convergent only in the strip < Re x < Q, otherwise it is an 
analytic continuation. In fact T(x) is an entire function of x with zeroes at a; = —nb — m/b 
and X = Q + nb + m/b with n and m non- negative integers. 

In the sense mentioned above the explicit results ( p..ll| ) and ( |1.20| ) constitute the exact 
construction of the Liouville field theory on a sphere. For example, the four-point function 
can be explicitly expressed in terms of the three-point function 



Gai 



a2,03)"4 



00 

i J C{ai,a2,Q/2 + tP)C{a3,a^,Q/2-tP)\F{A^^,A,x,)\^dP (1.22) 
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where the intergration is over the variety of physical states \vp) and F{Aa., A,Xi) is the 
four-point conformal block, determined completely by the conformal symmetry fl^^. In the 
four-point case, which we are considering now, the latter can be further reduced to a function 
of one variable, e.g., 

F(A„^,A,x,) = (1.23) 

(X4 - xi)-2^Ha;4 - 0:2)^^+^^-^^-^^ (X4 - X3)^^+^^-^^-^^(x3 - ^a)^^"^^-^^"^^ X 

\P\v) 

where 

(Xi - X2)(X3 - X4) 
(Xi - X4)(X3 - X2) 

Parameters are related to A^. as in eq. (|1.7|) and in the intermediate dimension A = 

gV4 + p'. 



2. The boundary Liouville problem 

The basic ideas of 2D conformal field theory with conformally invariant boundary were 
developed long ago mostly by J.Cardy pO| who also applied them successfully to rational 
CFT's, in particular to the minimal series ET], Here we'll try to apply these ideas to the 



Liouville CFT with boundary. 

A conformally invariant boundary condition in LET can be introduced through the fol- 
lowing boundary interaction 

v4bound = ^buik + j (|^<^ + fJ^Be'"^^ g^'^d^ (2.1) 
ar 

where the integration in ^ is along the boundary while K is the curvature of the boundary 
in the background geometry g. In what follows we consider only the geometry of a disk 
which can be represented as a simply connected domain F in the complex plane with a fiat 
background metric gab = Sab inside. The action is simplified as 

Abound = J (^(^a0)' + fie^'f^ d'x + J (^^0 + ^^e'^) (2.2) 
r ar 

where k is the curvature of the boundary in the complex plane. Typically the most convenient 
domain is either a unit circle or the upper half-plane. In the last case the boundary dT is 



3c 



Strictly speaking, (1.22) is literally correct only if Rea^ are sufficiently close to Q/2. Otherwise additional 



discrete terms in the r.h.s of (1.22) can appear due to certain poles of C bracking through the integration 
contour, see [|l9|l. 



6 



the real axis and one can omit the term hnear in in the boundary action ( |2.2| ). The price 
is again a "background charge at infinity", i.e., the same boundary condition on the field 
at infinity in the upper half plane 

(l){z,z) = -Q\og{zz) + 0{l) at |z| ^ oo (2.3) 

as in the case of the sphere. 

It seems natural to call the additional parameter /i^ the boundary cosmological constant. 
We see that in fact there is a one-parameter family of conformally invariant boundary condi- 
tions characterized by different values of the boundary cosmological constant /i^- Contrary 
to the pure bulk situation where the cosmological constant enters only as a scale parameter, 
the observables in the boundary case actually depend on the scale invariant ratio yu/yU^. 
For example, a disk correlation function with the bulk operators , ■ ■ ■ and the 
boundary operators (see below) , Bf^^ ■ ■ ■ B^^ scales as follows 

g{a,, ■ ■ - a., A, ■■■Pm)- (2.4) 

where F is some scaling function and we indicate only the dependence on the scale parameters 
H and /isQ. Our present purpose is to study this dependence. 

In the boundary case we have to introduce the boundary operators. In LFT the basic 
boundary primaries are again the exponential in boundary fields Bj^ = exp /?0. Their 
dimensions are 

Ap = P{Q-P) (2.5) 

To avoid any confusions we shall always use parameter a for the bulk exponentials and 
parameter jS in relation with the boundary operators. In general a boundary operator is not 
characterized completely by its dimension, because the conformal boundary conditions at 
both sides of the location of the boundary operator may be in general different. One has 
to specify which boundary condition it joins. Therefore in general we are talking about a 
juxtaposition boundary operator between, in our case, two boundary conditions with the 
parameters /i^i = yUi and yUsa = and denote it B^^^'^^x). 

To define completely the boundary LFT on the disk, i.e, to be able to construct an 
arbitrary multipoint correlation function including bulk and boundary operators, we have to 
reveal few more basic objects in addition to the bulk two- and three-point functions ( p..l8|) 
and ( |1.2CI| ) we already have. 



1. First is the bulk one-point function (we imply almost constantly the upper half-plane 
geometry) 

iV.(^)) = (2.6) 

\z — z\ 

In fig.|l|a it is drawn however as the one-point function in the unit disk. 

"'in the presence of boundary operators it is possible to impose different boundary conditions at different 
pieces of the boundary, each being characterised by its own value of /is. In this case the scaling function in 
(^.4|) may depend on several invariant ratios, see below. 
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(a) 



(b) 




(c) 




(d) 



Figure 1: 



2. Second, one needs the boundary two-point function 

/i2) 



\2Ar-, 



(2.7) 



which in general depends on two boundary cosmological constants /iiand ^2 (see fig.^b). 

3. The bulk-boundary structure constant, which determines the fusion of a bulk operator 
Va with the boundary resulting in the boundary operator B^^'^^. This is basically the same 
as the bulk-boundary two-point function (fig.|l]c) 



{V.{z)Bp{x)) 



R{a,P\fiB) 



\z — z\ 



Z — X\ 



\2Ar- 



(2.8) 



In fact the one-point function ( p.6|) is a particular case of this quantity with /3 = so that 
its introduction, however convenient, is redundant. 
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4. Finally, there is a boundary three-point function 



{^p, (xiji^^^ {x2}Bp.^ {X3)) - A1+A2-A3, 1A2+A3-A1, 1A3+A1-A2 l^-y) 

F12I F23I F3l| 

which in fact depends now on three different boundary parameters, /xi, fi2 and /is related to 
the corresponding sides of the triangle as shown in fig.|I|d. 

These three basic boundary structure constants, together with the bulk structure con- 
stants, allow in principle to write down an intermediate state expansions for any multipoint 
function. An instructive example is the bulk two-point function 

Ga,aM,Z2) = (V;,(^1)V;,(Z2)) (2.10) 

Joining these two operators together with the bulk structure constant we can reduce this 
quantity to the one-point bulk function and write down the following expansion 

- |2Ai-2A 



= ''^ ,4A, —C{a,,a,,Q/2 + iP)U{Q/2-iP)T[ |P|J 

\Zl - Z2\ J \ 012 012 J 



(2.11) 



where 



= (2,12) 

[Zl - Z2){Zi - Z2) 

is the projective invariant of the four points zi, Z2, zi and Z2 and T is the same four-point 
conformal block as enters the expansion of the four-point bulk function, see ( |1.22|jl.2B| ). 
Notice that while in that case it entered in a sesquilinear combination ( |1.22| ), here it appears 



linearly (J.Cardy |^). Expansion (p.ll|) is suitable appropriate if the bulk operators are 



close to each other, i.e., 77 0. Another representation is suitable the limit rj ^ 1 where 
the points zi and Z2 approach boundary and the bulk operators can be expanded in the 
boundary ones. This gives 



I - |2Ai-2A2 p ip / \ 

G'aia2 = ^^^^^^A^ / — ^(«i,Q/2 + ^m«2,g/2-^P)^ \p\i-v] 
\zi - Z2\ J \ ai a2 J 

(2.13) 

Equating these two expressions we see that the basic boundary quantities also must satisfy 
some bootstrap relations analogous to that in the bulk case. It is interesting to note, that 
there is another application of this relation. The conformal block itself, although being 
completely determined by the conformal symmetry, is it fact a complicated function which 
is not in general known explicitly. On the other hand it is important since it explicitely 



enters the conformal bootstrap equations [T^. Besides, one might expect that it encodes 
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some information about the structure of the representations of the conformal symmetry. In 
particular conformal block must satisfy the following cross-relation 



^f"^ \P\v] = f —K ( |P,P'Vf"' 1^11-^1 (2.14) 

\ a2 ' ' ' J J 2-K \ a2 ' J \ as 0:4 ' ' 'J 

with some cross-matrix K which determines the monodromy properties of the conformal 
block. Suppose now we've managed to find the basic quantities of the boundary Liouville 
problem, in particular the one-point function U (a) and the bulk-boundary structure constant 
R{a,(3). Then the crossing relation becomes a linear equation for the cross-matrix of the 
symmetric (i.e., 0:3 = ai and 04 = 02) conformal block, from where this matrix can be 
figured outQ. 



2.1. Bulk one-point function 

We start with the calculation of the bulk one point function U{a\fiB)- For this we apply 
the degenerate operator insertion, like above for the bulk two-point function. Consider the 
auxiliary bulk two-point function with the additional degenerate bulk field V_h/2{z) 

Ga,-b/2ix, z) = {V^{x)V^„2{z)) (2.15) 

Apply first the OPE aX z ^ x where the degenerate operator K_b/2 generates only two 
primary fields so that 

Ga,^b/2 = C+ {a)U{a - b/2)g+{x, z) + C_(a) U{a + b/2)g_{x, z) (2.16) 

where C±{a) are the special structure constants as given by the screening integrals and 
Q±{x^ z) are expressed through the special conformal blocks J-'±{x, z) related to these special 
values of parameters 

Q±{^.^) = ' , ' _,4A. ^±iv) (2.17) 



X 



where A_b/2 = -1/2 - Sb^/A and 

n = 

{z — x){z — x) 



V=['-fr''\ (2.18) 



In fact satisfies the second order differential equation. Therefore these special conformal 
blocks are solution to a second order linear differential equation and can be expressed in terms 



^In a recent paper an explicit expression for this matrix has been proposed on the basis of completely 
different approach 
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of the hypergeometric functions 

_ . , _/ -6/2 -6/2 I b 

r/"*(l - ri)-^'/^F{2ab - 1 - 26^ -b^, 2ab - b\ 7]) (2.19) 

_ , . _/ -6/2 -6/2 I ^6 
\ a a 2 
^i+b^-ab^-^ - r])-'''/^F{-b^, 1 - 2a6, 2 + 6^ - 2«6, r]) 

This is a particular case of more general conformal block with a degenerate operator 

l«2-^|r/) =v''''{l-Vr'F{AB,C,r]) (2.20) 





f -6/2 


"3 






0:4 




f -&/2 


"3 








where 







l«2 + 2 1^ ) = V '(1 - r7)''°3F(A - C + 1, 5 - C + 1, 2 - r^) 



A = -1 + 6(^2 + «3 + «4 - 36/2) 

5 = 6(^2 + tt3 - tt4 - 6/2) (2.21) 
C = 6(2^2 - 6) 

Now, as both operators approach the boundary, they are expanded in the boundary 
operators. It turns out that the degenerate bulk operator V_b/2 near the boundary gives 
rise to only two primary boundary families Bq and B_fy. The simplest thing is to find the 
contribution of Bq = I. The fusion of Va to the unity boundary operator is described by the 
quantity R{a, 0) = v{a) while the fusion of the field V-b/2 (y~b/2 boundary) is described by 
a special bulk-boundary structure constant R{—b/2, Q). It can be computed as the following 
boundary screening integral with one insertion of the boundary interaction —^b J Bh{x)dx 

R{-b/2,Q) = -22^^>B J {V^b/2{t)Bb{x)BQ{^))dx (2.22) 

27r/iBr(-l-262) 
- r2(-62) 

Comparing this with the behavior predicted by the bulk expansion ( |2.16| ) we find the follow- 
ing functional equation for the one-point function 

.U(a) = - 6/2) - ^-^ipl,Z^^U(a + i/2) (2-23) 



r(-62) ' ' r(-l- 262 + 26a) ' ' ' -f{-b^)T{2ba) 

(in the last term we used the bulk special structure constant C_(a) from (|1.16|) ). Equation 
is solved by the following simple expression 

U(a) = \ (nfi^(b^)) Y(2ba - b^)T f ^ - 1 - 1^ cosh(2a - Q)tts (2.24) 

6 \ b ¥ J 
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where the parameter s is related to the scale invariant ratio of the cosmological constants 

cosh^ nbs = ^ sin vrfo^ (2.25) 
/i 

Also this expression satisfies the dual functional equation provided the dual bulk cosmological 
constant ft is related to fi as before in ( |1.19|) while the parameter s is self-dual, i.e. the dual 
boundary cosmological constant ftB is defined as follows 

cosh^ = ^ sin — (2.26) 

It is remarkable enough that the expression ( p.24| ) automatically satisfies the "refiection 
relation" JTsf for the operator 

U{a) = U{Q ~a)D{a) (2.27) 
with the bulk Liouville two-point function (|1.18| ). If a corresponds to a physical state, i.e.. 



a = Q/2 + iP with P real, expression ( |2.24| ) reads 



U{P) = (vr^^(62))-^/^r(l + 2z6P)r(l + 2zP/6)^^^i^^ (2.28) 
^ iP 

This quantity is interpreted as the matrix element between a primary physical state \vp) 



from (|1.8|) and the boundary state {Bs\ created by the boundary action ( p.2|) 

f/(P) = {B,\vp) (2.29) 
It is natural that this matrix element satisfies the refiection relation 

U{P) = S{P)U{-P) (2.30) 
with the Liouville refiection amplitude ||13| 

O /OA f ,^2AX-2iP/fer(l+2iP/6)r(l + 2iP6) 

s,iP) = - {n,,ib^)) r(i-2.p/6)r(i-2.P6) ^'-''^ 

Of course the functional relation does not fix the overall constant so that it can be 
multiplied by any (self-dual) factor Uo{b) = Uo{l/b). In ( p.24| ) this factor is chosen in the 



way that all the residues in the "on-mass-shell" poles at 2a; = Q — nb, with n = 1, 2, 3, . . . 
are equal precisely to the corresponding perturbative integrals appearing in expansions in fj, 
and fiB, i.e. 

- V (-^^n-^^Br-'' X (2 32) 



k=0 



lmzi>0 



oo 

'Zk j dXi... dXn-2k {VQ^nb{'i)Vb{Zi) . . . Vb{Zk)BhiXi) . . . Pfe(x„_2fc))c 
— oo 
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where (■ ■ ■ )o is the correlation function w.r.t the upper half plane free field with fi = fiB = 
i.e., with free boundary conditions. Explicitely 



n 

nl - \-'2a^ TT I 

i=l i,j 



i>j i>j 



(2.33) 



In particular, the pure boundary perturbations in /x^ reproduce the Dyson integrals over a 
unit circle 



n 

.2 



r(l - nV) (2.34) 



Several remarks are in order in connection with the expression (|2.24 ) presented. 



1. Semiclassical tests. Consider the limit 6^0 while P in eq.( |2.28|) is of order h 
and s is of order h^^. In this limit the minisuperspace approximation is expected to work. 
Take the geometry of semi-infinite cylinder of circumference 2tt and consider the states on 
the circle. In the minisuperspace approximation one takes into account the dynamics of the 
zero mode 



1 



27r 



^0 = — (f){x)dx (2.35) 
27r Jo 



neglecting completely all the oscillator modes of field 0(x). The primary state \vp) is repre- 
sented now by the wave function 



M<Po) = w o-o/.^ K^^p/t (2./^e''^A (2.36) 



2 MbT'"'/' /:^.6,o 

T{-2iP/h) 

{Ky{z) is the modified Bessel function) which satisfies the minisuperspace Shrodinger equa- 
tion 

and has the following asymptotic at 0o ^ ~oo 

^p(0o) = e2*^<^» + 5f )(P)e-2*^<^» (2.38) 

where S^^\P) is the classical limit of the Liouville reflection amplitude (|2.31| ). Also it meets 
the normalization (|1.9|) 

roo 

#,(0o)V^p(0o) = 7r(5(P - P') (2.39) 
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In the approximation under consideration the boundary state {Bs\ wave function is simply 
related to the boundary Lagrangian 

^B.(0o) = exp(-27r/iBe^<^«) (2.40) 

The matrix element {Bs\vp) can be carried out explicitely 



oo 



(0o)V^p(0o)c^0o = T (vr/i/6') r {2tP/b) cos(27rPs) (2.41) 



CO 



and agrees precisely with the corresponding limit of ( p.28| ). Note, that this calculation is 



sensitive to the prefactor in eq.( |2.24| ) and confirms our choice Uo{b) = 1, at least in the limit 

2. Boundary length distribution. From the point of view of 2D gravity one can 
interprete the quantity 

/ = / expb(f) (2.42) 
Jar 

as the lenth of the boundary. Let Wa{l) be the boundary length distribution for the fluctuat- 
ing disk with the bulk cosmological constant fi and an insertion of the operator somewhere 
inside the disk. Then 

U{a\iiB) = / ^e-^^'W^{l) (2.43) 







Form the result (|2.24|) one finds explicitely 



where 



= (2.45) 

Sm TTO^ 



Compare ( p.44| ) with the minisuperspace distribution ( p.36| ). This result implies that the 



Shrodinger equation ( p.37|) in the logarithm of the scale (po = b ^log(//2'7r) (which is called 



sometimes the Wheeler-deWitt equation) does not hold only in semiclassical limit but is 
in fact exact with a suitable renormalizations of constants (see in this relation the paper 
where this equation first appeared in the context of the Liouville field theory, and also 



| 25| , p6l where similar expressions are obtained in the framework of random surface models). 



Let us present also the double distribution in the length ( p.42| ) and area A defined as 

HA 

— e-^^Z,(AO (2.46) 
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It is given by a rather simple expression 

ZM. I) = exp --r— -7, (2.47) 



hV{l + l/h^ -2a/h) \ 2A J 4Asin7r62 

3. "Heavy" a semiclassics. Consider again the hmit 6^0 but with large value of 
a = ri/b not nessesserily close to Q/2. Exact expression (|2.47|) gives in this limit 

where C is the Euler's constant and 

S,{A, /) = ^ + (1 - 2r]) (log i2A/l) + log(l - 2r]) - 1) (2.49) 

On the other hand the corresponding classical solution with the area A and the boundary 
lenght / reads for the classical field ip = 260 (we imply here the geometry of the disk < 1 
with the unit circle as the boundary) 

^^(.) ^ P{a-'-af 

[27i{a-^{zzyn -a{zzy-^)f 



where a is related to the follows 



AnA 



a' = l-^{l-2v) (2.51) 

(we imply here that rj < 1/2 and P > 47rA(l — 27]) so that a real classical solution exists). 
The classical Liouville action for this solution is readily carried out 

Sci{A, /) = ^ + (1 - 2r/) (log {2A/1) + log(l - 2r/) - 1) (2.52) 

and coincides with ( |2.49| ). In principle it might be possible to check the prefactor in ( ^.48 ) 
performing the one-loop correction. This is not yet done however. 

4. "Light" a semiclassics. Direct semiclassical calculation of the one-point function 
( p.24| ) is possible also in the case a = ba, with 6^0 and a fixed. In particular, one can 
calculate the semiclassical approximation to the function (|2.47|) by taking the saddle-point 
contribution to the corresponding functional integral over (p with fixed area A and boundary 
length I. In the present case a ^ b the exponential insertion does not affect the saddle-point 
configurations. The nature of these classical solutions depends on the relative value of A and 
P. Here we consider explicitely only the negative-curvature situation Arc A < P, in which 
case the classical configurations form an orbit under the action of SL{2,R). To be specific 
we adopt the upper half-plane geometry with the boundary at the real axis. Then generic 
classical solution (pc is obtained from the "standard" solution 

e'^^' (z, .-) dz d-z = ^f '^"^"' (2.53) 

{{z + i){z — I) — v^Y 
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by SL{2, R) transformation 

cz + d \c d 

here 



G'*z = ^^, G= eSL{2,R) (2.54) 



21A WP-AtiA 
P=p-2^^ P-2nA 

The semi classical approximation to the expectation value ( 2.47|) is then evaluated as an 
integral over this manifold of classical configurations, i.e. 

(e^ba<t>^^^^ = Ne-^"^^' J dfi{G) {e^'"^^{z,z)Y (2.56) 

where is the classical action ( |2.52D with r] = 0, d^{G) stands for the SL{2, R) invariant 
integration mesure and the factor combines the determinant of zero modes and the con- 
tributions of positive modes to the gaussian integral around given classical solution. It is 
important to note that while N can very well depend on A/P, it carries no dependence on 
a, i.e. all the a dependence of the one-point function in this approximation comes from the 
integral in ( p.56|) . 

The integrand in (|2.56|) can be simplified by a shift of the integration variable G ^ GzG, 
where Gz is any fixed (2;-dependent) SL{2, R) transformation which maps the point z to the 
point i in the upper half-plane; this gives for the integral in ( |2.56 ) 

(^-^Y j ^-"(^^ + 6' + + + 2) - z/2(c2 + t;2))-2._ (2.57) 

To evaluate this integral one can introduce the following coordinates on the group manifold 
of SL{2,R), 

2yy-x{y + y) y + y-2x 1 
a/d = i ; bld = x\ cld = % — ; d 



y-y y-y v$^m Y\y-x\ 

(2.58) 

where x is real and y and y are complex conjugate with y > 0. The invariant mesure 
takes the form 

, , ^, 2i dx d'^y 

MG) = 7 12 2.59 

[y-y) \y-x\^ 

and the integral in ( |2.57| ) can be written as 

y-yy^'' {{y + ^){y-^)-ur''P^ (2.60) 



2i J \y — x\ 
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This integral is readily evaluated and one obtains for ( p.56| ) 

N e^i^^^' 2.61 

2a-l\z~z I J ^ ' 

where the factor N = ttP N/A does not depend on cr; as is mentioned above its determination 
requires analysis of the fluctuations around the classical configurations which we did not 
perform. The cr-dependent part in ( p.61|) agrees with 6 — >■ limit of (|2.47|) . 



5. Boundary state. Once the function U{P) is constructed, the boundary state {Bs\ 
can be written down explicitely 

r°° dP 

(5.1= / —U{P){P\ (2.62) 



where the so called Ishibashi states ||2 



are designed in the way to match the conformal invariance of the boundary. Since the 
combination U{P) (P| is invariant w.r.t. the reflection P —>■ —P one can extend formally 
the integral ( |2.62|) to the negative values of P and write 

roo ip 

{Bs\= e2-^^w(P)(P| (2.64) 

J-oo ^TT 

where 

uiP) = ^ (7r/i7(62)) r (1 + 2^6P) r (1 + 2iP/b) (2.65) 
iP 

It is natural to call u{P) the boundary state wave function. Note that the state (P|, although 
consistent with the conformal invariance, does not correspond to any conformal boundary 
state, i.e., to a state created by a local conformally invariant boundary condition. However, 
it can be constructed as a linear combination of boundary states. In view of eq.j^M) we 
can write down 

^ (P| = r (P,| ds (2.66) 

This equation allows to single out a conformally invariant state containing only one primary 
state {vp\ and its descendents. In finite dimensional situation of rational conformal field 
theories this trick has been friquently used by J.Cardy [Q. 
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3. Boundary two-point function 



In this section the boundary two-point function d{(3\fii, ^2) of ( p.7|) will be derived. To this 



purpose we apply basically the same Techner's trie which has been used in the first section 
to determine the bulk structure constants. Considering the boundary operators Bpi^x) we 
find that all the operators B_nb/2{x) (and also of course the dual fields B_n/2b) with n = 
0, 1, . . .are degenerate, i.e., count primary states among their descendents. A complication 
here is that not all of these "null vectors" nessesserily vanish, contrary to what happens in the 
bulk situation. For example simplest non-trivial degenerate boundary operator B'^^^J^ (form 

now on we shall denote the exponential boundary operators -B^^*^ = ef^^^^ instead of B^^^'^ 
having in mind the relation ( |2.25| )) in general does not satisfy the second order differential 
equation. This means that the null-vector in the corresponding Virasoro representation is 
some non-vanishing primary field and therefore the second order differential equation has 
some non-zero terms in the right hand side. This effect can be already seen at the classical 
level where the upper half-plane boundary Liouville problem is reduced to the classical 
Liouville equation for the field = 



dd^ = 27r/i62e^ (3.1) 
in the upper half-plane with the boundary condition 

i{d - d)if = AiTfiBb^e^/^ (3.2) 
at the real axis. The boundary value of the classical stress tensor can be easily computed 

^ci(a;) = -j^^l + \^xx + nb^Tifilb'^ - /i)e'^ (3.3) 

The boundary operator B^_^^^^ in the classical limit reduces to the boundary value of exp(— (^/4) 
for which we have 



3<p/4 



(3.4) 



In the right-hand side there is a primary Virasoro operator exp{3(p/4)s^s which has exactly 
the same dimension as the null-vector in the corresponding degenerate representation. It 
is interesting to note that there is a unique relation between the cosmological constants 
7r/i^6^//i = 1 where the r.h.s vanishes and this operator satisfies homogeneous linear dif- 
ferential equation. This effect holds on the quantum level too: if the boundary and bulk 
cosmological constants are related as 

l = ^tan^ (3.5) 

the second order differential equation holds for the boundary operator B^^^i^ (see the remark 
in the concluding section). 
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Here we are interested in the general situation where this operator is of no use since it 
does not always satisfy the second order differential equation. It happens however that the 
next degenerate boundary operator B'^l do satisfy the third-order differential equation when 
placed between identical boundary conditions. Therefore it can be used in our calculations 
instead of -B_b/2- As in the bulk, the differential equation predicts the following truncated 
OPE of this operator with any exponential boundary primary 

B'4Bf = c4Bp_,] + co[Bp] + c_[fis+5] (3.6) 

where Ca{P) are again the special boundary structure constants, which can be calculated as 
certain screening integrals. Considering again the auxilary three-point boundary function 
with a B_h insertion one figures out immediately that 

d(0 + l\,u,,) ^ (3.7) 

d{fJ\ni,H2) 

The structure constant c_(/3) can be evaluated as a combination of scrining integrals. 
These are of two tipes: a volume screening by the bulk Liouville interaction term e^'"^ 



(v) 

C_ 



lmz>0 

/ d'^z- 



lmz>0 

J2, 



|1 - z\ 
, 1 46/3 I -1 262 

J \z\ \z — z\ 

lmz>0 

and two boundary screenings e^^ related to the boundary interaction 



I I rfxidx2(e^*(-)e^<^(-)i?--(0)i?T(l)i??-}3-6(oo)> (3.9) 



/^i/^J- f f A^. A^. I(l-^l)(l-^2)| 



/ / dxidx2 



2b' 



,2 



oil "'•^l"-^^ 262 26/3 

^ -'Ci JCj \Xi - X2 0:1X2 



where the contours Ci are chosen as in fig.|] while /ij are the corresponding values of the 
boundary cosmo logical constant, as it is also indicated in in the same figure. Both contribu- 
tions can be carried out explicitely and we have 

c_(/?) = cL"^ + cL^) 

= [-/isin2(27r6/3) + sin7r62 (/^^ + ^2 _ 3^^^^ cos(27r6/3))] Iq{(3) (3.10) 
= AjjL sm no sm vro- 



2 / V 2 
sm TTO sm TTO io\p) 
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where 



7(1 + 



TT 



-r(l - 2bp)T{2bp - l)r(l -b^ - 2b/3)T{2b/3 - b^ - 1] 



and Si and S2 are again related to /ii and ^2 as in eg. (|2.25|) . 



(3.11) 



(1) 







^2 



^2 



Figure 2: Contours of integration Ci in (|3.9| ) together with the corresponding values of the 
boundary cosmological constants /i,. 



To construct a solution to the functional equation { \i.T\ ) with ( |3.10| ) we need more special 
functions. First one is what is sometimes called the g-gamma function. Here we denote it 
S(x). It is self dual with repect to b ^ 1/b and satisfies the following shift relations 



S(x + b) = 2 sin(7r6x)S(x) 
S{x + l/b) = 2sin(7ra;/6)S(a;) 



(3.12) 



It has zeroes a.t x = Q + nb + 'm/b and poles at a; = —nb — m/b ( m and n are non- negative 
integer numbers). In the strip < Re a; < Q the following integral representation is allowed 



log S(x) 



dt 
T 



sinh(Q-2x)t {Q/2-x) 



2sinh(6t) sinh(t/6) 



(3.13) 



With this definition it satisfies also the "unitarity" relation 

S{x)S{Q-x) = 1 



(3.14) 



It is also convenient to introduce a self-dual entire function G(x) which contains only zeroes 
at a; = —nb — m/b, m,n = 0,1,2, .. . and enjoes the following shift relations 



hl/2-bx 

G{x + b) = -^^T{bx)G{x) 
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Ux/b-1/2 

G(x + 1/b) = ^ r(x/&)G(x) 



(3.15) 



27r 
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This function is "elementary" in the sense that both T(x) from eq.( |1.21| ) and S(a;) are simply 
expressed in G{x) 



T(x) = G{x)G{Q-x) 
GiQ-x) 



Six) 



G{x) 



The integral representation which is valid for all < Re x reads 

oo 

Q-Qt/2_^-xt {Q/2-x^^ 



log G{x) 



dt 
T 



+ 



Q/2~x 



(l-e-^*)(l-e-*/^) ' 2 



With this function one can easily construct a solution to ( |3.7| ). 

G(Q-2/3)G^i(2/3-Q) 



(3.16) 



(3.17) 



d{/3\si, S2) 



„ , (Q-2/3)/2fe 

7r/i7(&')&'-'' 



S{(3 + t{s, + S2)/2)S{f3 - t{s, + S2)/2)S{(3 + t{s, - S2)/2)S{f3 - t{s, - S2)/2) 

(3.18) 



exp 



dt 
T 



sinh(Q — 2P)t cos sit cos S2t {Q — 2/3) 
sinhfet sinht/6 t 



This solution satisfies also the "dual-shift" relation analogous to ( |3.7| ) so that ( |3.18D is the 
unique self-dual solution to ( p.7|) . It is of course possible to express the ratio of two G- 
functions in terms of S-function times some ordinary F-functions. We prefere to present 
d{f3\si,S2) in the form ( |3.18| ) to make obvious the "unitarity" relation 



d{(5\si, S2)d{Q - (3\si, S2) = 1 

Note, that an overall independent of /5 constant which is allowed by 
completely fixed by ( |3.19D . 

4. Concluding remarks 



(3.19) 
and its dual is 



Eq. (|3.5|) together with the structure of singularities of the two-point function ( p.l8|) 
drop a hint at the suggestion that the level 2 degenerate boundary operator B'^_^^^^{x) 
has a vanishing null vector if and only if Si — S2 = ±ib or Si + S2 = ^ib (the second 
condition is requred by the symmetry of boundary conditions w.r.t. s — > — s). Let us 
verify this suggestion on the three-point function with one boundary field i?^'^^*''(x), 
i.e., consider 



BrHx,)B%J,^^\x)Blf^;'^^ix2) 



(4.1) 



21 



Under our suggestion _B^'^^*^(x) satisfies a second order differential equation in x and 
therefore lias special operator product expansion with any Bp 

(4-2) 

Then, exactly the same trick which led to eq. (|3.7|) gives the following shift relation 

d{P\s,s') 



TyS.sUb j-)S±ib,s' (±) I j-)S,s' 



(±)[rs,s' 

l3-b/2i 



d{(3 + b/2\s±ib, s') 



As usual we adopt the structure constant with no screenings requred c, 
structure constant CL^ is given by the integral 



(±) 



J — oo 



oo) ) dx 



This integral is evaluated quite easily (unlike (|378|) or (|3.9| )) 



r(l - 2bp)T{2b(3 - 62 - i)r(l + 62 

TT 



(4.3) 



1. The 



(4.4) 



(4.5) 



jji sin7r(62 — 26/3) + jj,^^^ sin7r(26/5) — sin rcb'^ 

where is determined by the relation 

cosh7r6(si ± ib) = /i^^'a/ simrb'^ / fi 
After some simple algebra we obtain 

\ 1/2 
7r7(— 62) J 

sin TTb {13 =F ib{si + S2) /2) sin 7r62 (/3 =F ib{si - S2) /2) 



r(l - 26/5)r(26/5 - 62 - 1) X 



(4.6) 



(4.7) 



It is easy to see that the two-point function ( |3.18| ) satisfies both relations ( |4.3| ). After 
this support one may suggest further that any degenerate field -Bl'^{,/2 vanishing 
null-vector (and therefore has truncated operator product expansions if s — s' = ibk or 
s + s' = ibk with k = —n/2, —n/2 + 1, —n/2 + 2, . . . , n/2, in close analogy with the 
fusion rules for degenerate bulk fields. 

The boundary two-point function (|3.18| ) is readily applied as the reflection coefficient in 
the reflection relations for the one-point function of an exponetial boundary operator in 
the boundary sin-Gordon model. The latter is defined by the following two-dimensional 
euclidean action 



A 



bsG 



d X 



An 



{da4>f - 2/i cos(2/30) 



-2^iB / cos(/3(0-0o)) (4.8) 



ar 
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where the bulk part of the action is integrated over a half-plane T so that the boundary 
dr is a strainght line. For the moment /3 denotes the standard sin-Gordon coupling 
constant. Apart from it the boundary model depends of three parameters /i, /xb and 
00 HH- The dimensional parameters fi and fiB can be given a precise meaning by spec- 
ifying the normalisation of the composite fields they couple to. As these operators are 
combinations of exponentials it suffices to specify a normalisation for the exponential 
fields in the volume and at the boundary. Here we adopt the conventional normali- 
sation of these fields (see e.g. [EH]) corresponding to the short distance asymptotics at 



\x 







e^''"''(x)e-^''"'' 



e""^^(x)e 



(y) 



1 



\x - y\ 
1 



l4a2 



\X 



y\ 



|2a2 



+ 



+ 



for the volume fields 



for the boundary ones 



(4.9) 



Here we present only the result for the one-point function of the boundary operator 
Gb{(i) = (exp(m0B)) which reads 



TT/i \2(1 



' 9o{a)gs{a)9A{a) 



(4.10) 



where 



log 5^0 (a) 

^og gs{a) 
log gA{a) 



dt 
T 



2sinh^(a/3t) (e 



,{l-/32)t/2 



cosh(t/2)cosh(/?2t/2) - 1 



sinhtsinh(/52t) sinh((l - (3^)t) 



— a e 



rftsinh^(a/?t) (2 - cos(2;zt) - cos{2z*t)) 
1 sinhtsinh(/?2t) sinh((l - p'^)t) 
dtsmh{2a/3t) (cos(2zt) - cos(2;z*t)) 



(4.11) 



t sinh t smh{(3H) cosh((l - (3^)t) 
where the complex number z is related to the parameters of the model 



as 



cosh nz 



sin nP 



(4.12) 



and z* is the complex conjugate to z. The details and some applications will be 
published elsewhere. 



In the next publication |^ we will present an explicit expression for the bulk-boundary 
structure constant (|2.8D. Equation ( ^2.66] ) then permits to resolve the system ( |2.11| ), 
( |2.13| ) and ( p.l4[ ) for the cross-matrix and obtain an explicit expression for a special 

Ci2 0-2 



case of symmeric cross-matrix K 



P,P' 
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The random lattice models of 2D quantum gravity allow in many cases to find ex- 
plicitely the partition functions of minimal models on fluctuating disk with some bulk 



and boundary operators inserted [25, 26]. Detailed comparison with the Liouville field 



thery predictions seems quite interesting. The work on that is in progress. 
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